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Ultrasound  Perfusion 


INTRODUCTION  (excerpts  from  the  original  proposal) 

Nature  of  the  Problem 

Mammography  has  proven  to  be  reliable  as  a  screening  tool  for  breast  cancer.  However, 
the  specificity  of  mammography  for  breast  cancer  may  be  as  low  as  10%  as  evidenced  by  the 
number  of  biopsies  recommended  compared  to  the  number  of  cancers  confirmed  (Moskowitz  and 
Gartside,  1982).  Therefore,  of  significant  health  care  benefit  would  be  a  method  used  in 
conjunction  with  mammography  which  could  reduce  the  number  of  biopsies  required  while 
maintaining  or  improving  survival  rates.  In  terms  of  health  care  costs,  Adler  et  al  (1990) 
estimated  that  more  than  half  of  the  mammographic  screening  costs  for  breast  cancer  are  the  result 
of  required  biopsies  or  excisions  due  to  low  specificity.  Again,  an  alternative  method  for 
establishing  the  type  of  lesion  present  would  be  of  great  benefit.  Any  improvement  in  early 
detection  of  breast  cancer  is  clearly  important  as  survival  is  significantly  improved  when  cancers 
are  detected  in  the  0.5  to  1  cm  diameter  size  range  (Axtel  et  al,  1976). 

Background  of  Previous  Work 

Doppler  ultrasound  has  been  investigated  and  found  some  utility  in  the  detection  of  cancer 
(Bums  et  al,  1982;  Minasian  and  Bamber,  1982;  Boyd  et  al,  1983)  .  Although  the  specificity  has 
been  high  for  cancer  recognition,  at  10  MHz  frequencies  used  in  these  studies,  the  sensitivity  of  the 
technique  has  not  and  the  use  of  other  Doppler  techniques  have  had  mixed  specificity  and 
sensitivity  results  ( Rubin  et  al,  1987;  Adler  et  al,  1988;  Jellins,  1988;  Cosgrove  et  al,  1990;  Adler 
et  al,  1990).  The  advent  of  stable  ultrasound  contrast  agents  and  the  development  of  Doppler 
power  mode  imaging  are  clearly  innovations  which  could  improve  the  sensitivity  of  ultrasound 
techniques  for  cancer. 

Abnormal  accumulation  of  contrast  during  the  arterial  phase  of  contrast  transit  and  a 
residual  concentration  of  contrast  agent  in  breast  cancers  has  been  observed  with  conventional 
angiography,  but  without  exceptionally  high  specificity  (Feldman  et  al,  1967;  Kaushik  et  al,  1975; 
Sakki,  1974),  with  x-ray  computed  tomography  (Chang  etal,  1982)  and  with  digital  subtraction 
angiography  (Flynn  et  al.,  1984).  Gadolinium  MRI  is  now  looking  quite  positive  as  a  means  of 
discriminating  likely  breast  cancers  and  even  detecting  mammographically  occult  malignancies 
(Harms,  et  al.,  1993,  Heywang-Kobrunner,  et  al.,  1992).  With  the  advent  of  ultrasound  contrast, 
such  contrast  studies  should  now  be  possible.  The  desired  IV  administration  of  the  agent  however 
results  in  dilution  of  the  contrast  and  the  loss  of  a  well  resolved  bolus.  If  such  problems  can  be 
overcome,  the  prospect  for  ultrasound  contrast  detection  of  breast  cancer  would  be  greatly 
improved  due  to  the  increased  signal  strength  from  the  smaller,  more  diagnostic  vessels.  Doppler 
power  imaging  has  also  demonstrated  an  improved  signal  to  noise  ratio  over  conventional  color 
Doppler  flow  iinaging.  Images  being  made  in  this  mode  seem  to  indicate  one  can  expect  such  a 
modality  would  improve  on  the  sensitivity  of  ultrasound  for  small  vessels. 

Finally,  recent  work  (Weidner  et  al.,  1991;  Horak,  1992)  concerning  angiogenesis  has 
placed  renewed  emphasis  on  the  blood  flow  in  the  region  of  cancerous  lesions.  The  potential  of 
ultrasound  techniques  described  in  this  proposal  to  measure  regional  perfusion  suggests  that 
information  can  be  derived  by  noninvasive  means  using  a  relatively  inexpensive  imaging  modality. 

Purpose  of  Present  Work 

The  present  research  is  design  to  develop  techniques  for  ultrasonic  measurement  of  tissue 
perfusion  in  the  breast.  In  combination  with  other  information  such  as  mammographic 
examination,  the  techniques  proposed  should  have  application  in  the  detection  of  cancerous  lesions 
or  the  monitoring  of  therapeutic  interventions  which  affect  the  vascular  flow  of  tumor  tissue.  The 
diagnostic  specificity  for  cancer  could  be  increased  with  these  techniques  leading  to  improved 
patient  management  while  using  an  imaging  modality  which  is  generally  less  expensive  than  many 
alternatives. 
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Methods  of  Approach 

This  study  examines  the  potential  of  three  related  Doppler  techniques  for  measuring  tissue 
perfusion.  1)  Acoustic  fields  can  be  used  to  modulate  a  continuous  infusion  of  bubble-based 
echocontrast  agent  to  estimate  wash  in-wash  out  time  courses  for  the  agent.  2)  The  decorrelation 
time  of  the  Doppler  power  signal  can  measure  the  regional  perfusion  based  on  the  rate  at  which 
blood  moves  out  of  the  Doppler  sample  volume.  3)  The  combination  of  Doppler  signal  power  and 
velocity  will  be  used  to  approximate  the  tissue  perfusion  where  the  Doppler  signal  power  is 
normalized  to  that  found  in  a  large  vessel  near  the  region  of  interest.  Development  of  these 
techniques  will  be  through  ex  vivo  and  in  vivo  experimentation  using  simple  and  physiologically 
relevant  flow  phantom  and  animal  studies.  The  latter  two  techniques  described  above  will  be 
examined  in  limited  clinical  trials  within  the  time  frame  of  this  study. 

DESCRIPTION  OF  RESEARCH  PROGRESS 

The  following  is  a  list  of  relevant  publications  which  have  resulted  from  this  work. 

Peer  Reviewed  Journal  Articles 

Adler  RS,  Rubin  JM,  Fowlkes  JB,  Carson  PL,  Pallister  JE:  Ultrasonic  Estimation  of  Tissue 
Perfusion:  A  Stochastic  Approach,  Ultrasound  Med.  Biol.  Vol.  21  (4),  493-500,  1995. 

Rubin,  JM,  Adler  RS,  Fowlkes  JB,  Spratt  S,  Pallister  JE,  Chen  JF,  Carson  PL:  Fractional  moving 
blood  volume  estimation  using  doppler  power  imaging.  Radiology,  accepted,  1995. 

Abstracts  and  Proceedings 

Fowlkes  JB,  Gardner  EA,  Carson  PL,  Ivey  JA,  Rubin  JM:  "Acoustic  Interruption  of  Ultrasound 
Contrast  Agents  for  Blood  Flow  Evaluation",  Procs.,  39th  Annual  Conv.  of  the  American  Institute 
of  Ultrasound  in  Medicine,  San  Francisco,  CA,  March  26-29,  1995. 

Rubin  JM,  Adler  RS,  Fowlkes  JB,  Spratt  S,  Pallister  JE,  Chen  J-F,  Carson  PL  (1995)  Fractional 
Moving  Blood  Volume  Estimation  Using  Doppler  Power  Imaging,  Twentieth  International 
Symposium  on  Ultrasonic  Imaging  and  Tissue  Characterization,  Arlington,  VA,  Program  and 
Abstracts,  June  7-9,  1995. 

Fowlkes  JB,  Gardner  EA,  Carson  PL,  Ivey  A  and  Rubin  JM:  Use  of  Contrast  Interruption  in  the 
Measurement  of  Blood  Flow.  Symposium  on  Ultrasonic  Imaging  and  Tissue  Characterization, 
Arlington,  VA,  June  7-9,  1995. 

Fowlkes  JB:  Progress  in  Perfusion  and  Moving  Fractional  Blood  Volume  Estimation  Using 
Doppler  Imaging.  Ninth  International  Congress  on  the  Ultrasonic  Examination  of  the  Breast , 
Indianapolis,  IN,  Sept.  28-October  1,  1995. 

Chen,  JF,  Fowlkes,  JB,  Carson,  PL  and  Moskalik:  Determination  of  Scanhead  Motion  Using 
Decorrelation  of  Speckle,  AAPM  37th  Annual  Meeting  and  Exhibition,  Boston,  MA,  July  23-27, 
1994. 


As  is  indicated  above,  this  research  project  is  divided  into  three  components.  Progress  in 
each  of  these  areas  will  now  be  discussed  in  turn.  Two  of  the  research  areas  have  been  combined 
in  this  writing.  This  method  of  relating  the  progress  information  is  chosen  for  presentation 
purposes  in  this  report  given  the  considerable  interaction  between  the  different  aspects  of  the 
project. 
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'  Contrast  Manipulation 

In  the  original  proposal  for  the  first  year’s  effort,  we  had  proposed  a  series  of  ex  vivo 
measurements  to  examine  the  parameters  which  control  contrast  manipulation  using  acoustic  fields. 
This  was  based  on  our  previous  observations  made  using  flow  in  tubes.  However,  we  were 
concerned  that  considerable  effort  could  be  spent  examining  the  process  under  an  initial  set  of 
conditions  which  might  not  be  sufficiently  like  those  found  in  vivo.  We  continue  to  see  some 
differences  in  the  case  contrast  agents  interrupted  in  blood  and  in  water  during  ex  vivo 
experiments.  Therefore,  we  conducted  two  initial  in  vivo  experiments  in  conjunction  with  our 
other  on  going  research  in  the  area  of  ultrasound  contrast  agents.  The  following  is  a  description  of 
those  experiments  and  the  results  of  one  of  these  is  presented. 

The  first  experiment  was  performed  in  a  rabbit  being  a  VX2  carcinoma  in  the  right  thigh. 
This  animal  was  being  used  to  evaluate  an  experimental  ultrasound  contrast  agent,  MRX-115 
(ImaRx  Pharmaceuticals,  Tucson,  AZ).  The  agent  is  a  perfluorocarbon  gas-filled  microbubble 

which  is  stabilized  with  a  lipid  shell.  It  has  a  mean  diameter  of  2.5  |im  and  a  bubble  concentration 

of  1.5  X  10^  bubbles/ml  with  99%  of  those  bubbles  being  <  10  pm  in  diameter.  Figure  1  is  a  size 
distribution  for  the  bubbles  as  measured  using  a  Accusizer  optical  sizing  system.  For  these 

experiments,  a  MRX-1 15  dose  of  10  pl/kg  was  used  for  each  injection  into  the  ear  vein  through  a 
catheter  and  followed  by  a  saline  flush.  The  animal  was  anesthetized  using  Ketamine  (5  mg/kg) 
and  Rompun  (35  mg^g)  and  the  right  leg  shaved  over  the  inner  (over  the  femoral  artery)  and  outer 
(over  the  VX2  carcinoma)  thigh.  Depilatory  creme  was  applied  to  the  skin  surface  to  remove 
residual  hair  and  the  skin  was  washed  with  a  surgical  scrub  and  rinsed  with  a  1:200  dilute  solution 
of  PhotoFlo  (Eastman  Kodak,  Rochester,  NY),  a  surfactant  to  act  as  a  wetting  agent.  The  animal 
was  then  placed  under  the  ultrasound  system  shown  in  Fig.  2  with  the  coupling  bath  over  the 
inner  thigh.  This  system  was  developed  for  transcutaneous  bubble  generation  (Ivey  et  al,  1995) 
and  uses  a  high  powered  coaxial  transducer  consisting  of  a  piezoceramic,  67  mm-diameter  outer 
element  used  for  producing  high  amplitude  fields  and  a  10  mm-diameter  inner  element  of  the  same 
material  used  for  pulse-echo  alignment.  The  outer  element  operated  at  a  center  frequency  of  1.8 
MHz  has  a  25%  bandwidth  and  a  focal  length  of  6  cm.  The  transducer  focal  beam  width  and  depth 
of  focus  were  1 . 1  mm  and  7.8  mm,  respectively  at  the  -6  dB  level.  A  single,  sinusoidal  tone  burst 
was  employed  with  durations  up  to  250  ms.  The  transducer  was  suspended  vertically  in  a  water 
standoff  tank  filled  with  degassed  water,  filtered  at  0.2  pm.  A  0.1  mm  thick  polyethylene  sheet 
was  used  for  the  acoustically  transparent  bottom  window  of  the  water  standoff. 

The  focal  point  of  the  power  transducer  was  aligned  at  approximately  the  center  of  the 
femoral  artery  lumen  using  a  Diasonics  Spectra  VST  (Diasonics  Ultrasound,  Milipitas,  CA)  and  its 
7.5  MHz  sector  format  endocavitary  probe  which  was  rigidly  fixed  to  the  transducer  so  that  the 
axis  of  the  power  transducer  was  in  the  image  plane.  The  image  of  a  mechanical  pointer  placed  on 
the  power  transducer  face  was  used  to  mark  the  location  in  the  image  plane  corresponding  to  the 
focus.  The  femoral  artery  was  identified  with  the  diagnostic  scanner  in  Doppler  mode  and  the 
combined  transducer  assembly  was  then  positioned  so  that  the  femoral  artery  was  at  the  point  in  the 
image  marking  the  focal  point  of  the  power  transducer. 

An  HP  3314  Function  Generator  (Hewlett  Packard,  Everette,  WA)  was  trigger  manually  to 
produce  the  gating  pulse  to  control  the  length  of  the  RF  burst  output  by  the  Wavetek  Model  50 
function  generator  (Wavetek  Inc.,  San  Diego,  CA).  The  RF  burst  was  amplified  by  an  ENI A300 
RF  amplifier  (ENI,  Inc.,  Rochester,  NY)  and  connected  to  the  1.8  MHz  transducer  through  a 
matching  circuit.  A  separate  10  MHz  linear  array  on  the  Diasonics  scanner,  operating  typically  at  6 
MHz  for  Doppler  modes  was  placed  over  a  distal  artery  in  the  VX2  carcinoma  and  used  in  the 
spectral  Doppler  mode. 
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Figure  1  -  Size  distribution  for  the  MRX-1 15  contrast  agent  used  for  in  vivo  contrast  interruption. 


Figure  2  -  Transcutaneous  contrast  interruption  setup  using  high  power  transducer  to  halt  the  flow 
of  contrast  injected  IV  in  the  ear  vein  of  a  rabbit.  The  targeting  transducer  was  the  Diasonics 
endocavitary  probe. 
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Figure  4  -  Example  of  transcutaneous  contrast  interruption  in  vivo  using  80  cycle  (~40  p,sec)  burst. 
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<  A  series  of  contrast  injection  were  performed  while  varying  amplitude  and  burst  length  of 
the  contrast  interruption  field.  Figure  3  shows  an  example  of  the  result  for  a  20  msec  burst  at 
approximately  20  MPa.  After  the  second  peak  systolic  portion  of  cardiac  waveform,  a  large 
vertical  band  is  seen.  This  is  an  interference  artifact  from  the  firing  of  the  high  power  transducer. 
During  the  third  cardiac  cycle,  an  abrupt  decrease  in  the  brightness  of  the  Doppler  spectral  display 
is  seen.  The  amplitude  of  the  display  had  been  set  so  as  not  to  saturate  in  the  presence  of  the 
ultrasound  contrast  agent.  Therefore  the  decrease  that  is  seen  is  down  to  the  level  of  the  pre¬ 
contrast  signal  coming  from  blood.  Note  that  the  decrease  comes  after  the  interference  effect 
owing  to  the  transit  time  for  the  negative  bolus  produced  in  the  femoral  artery  upstream.  What  was 

extraordinary  was  that  burst  length  could  be  reduced  to  as  low  as  40  cycle  (~20  |xsec)  and  the 

effect  could  still  be  seen  from  a  single  application.  Figure  4  is  the  result  for  an  80  cycle  (~40  jxsec) 
burst.  Note  that  this  time  the  interference  from  the  firing  of  the  high  power  transducer  is  not  seen 
because  it  was  so  short  on  this  time  scale.  However,  the  decrease  in  Doppler  signal  strength  is 
present  again  lasting  approximately  1  cardiac  cycle.  This  is  consistent  with  the  possibility  that  the 
shorter  pulse  is  sufficient  to  create  the  negative  bolus  in  the  femoral  artery  and  this  bolus  is  swept 
downstream  with  the  next  cardiac  cycle.  The  cardiac  cycle  duration  would  perhaps  ultimately 
control  the  duration  of  the  negative  bolus. 

As  part  of  another  project  exploring  the  use  of  high  power  ultrasound  to  transcutaneously 
generate  bubble  in  vivo  ,  we  were  also  able  to  perform  a  single  canine  experiment  to  further 
investigate  contrast  interruption.  The  experimental  protocol  was  essentially  the  same  as  that  above 
with  the  following  exceptions.  The  acute  animal  was  initially  anesthetized  with  suritol  (3mg/kg), 
intubated  and  anesthesia  maintained  by  ventilation  with  3%  enflurane.  The  contrast  interruption 
was  monitored  in  the  femoral  arteiy  immediately  downstream  from  the  interruption  site.  As  the 
animal  had  no  carcinoma  in  the  thigh,  the  site  selected  was  based  on  the  constraints  of  the  other 
experiments  being  performed.  Because  of  the  proximity  of  the  contrast  interruption  to  the 
observation  site  it  was  difficult  to  detect  the  negative  contrast  bolus,  however,  there  was  some 
indication  from  observations  during  the  experiment  that  contrast  flow  was  being  disrupted. 
Further  examination  of  the  data  has  not  yet  been  performed  to  confirm  these  observations  but 
should  be  available  at  a  later  writing. 

In  addition  to  these  specific  advances  in  contrast  interruption,  research  arrangements  with 
Diasonics  Ultrasound  have  resulted  in  some  techniques  which  we  will  be  able  to  apply  in  the 
coming  year.  Harmonic  imaging  is  now  available  on  the  Diasonics  Spectra  which  will  should 
provide  improved  signal  to  clutter  noise  ratios  for  this  research.  Harmonic  imaging  uses  the  B- 
mode  or  Doppler  imaging  modes  of  the  scanner  to  transmit  at  one  frequency  and  receive  at 
harmonics  of  that  frequency.  The  response  of  bubbles  at  the  second  harmonic  of  the  driving 
frequency  provides  signal  strength  which  is  orders  of  magnitude  greater  than  that  coming  from 
tissue.  Through  this  project  and  a  subcontract  to  an  NIH  SBIR  subcontract  we  have  measured 
second  harmonic  emissions  from  a  number  of  contrast  agents,  including  MRX-1 15.  Figure  5  is  a 
power  mode  image  of  a  canine  kidney  in  vivo  before  and  after  injection  of  the  MRX-1 15.  The 
dramatic  signal  increase  in  the  imaging  mode  as  the  agent  passes  through  indicates  that  this  agent 
should  provide  excellent  interruption  capabilities  for  regional  perfusion  estimates.  Experiments  are 
now  underway  to  test  the  enhancement  seen  in  the  second  harmonic  mode.  Our  initial 
measurements  of  the  second  harmonic  emission  indicate  that  the  agent  should  provide  reasonable 
second  harmonic  emission  as  seen  in  Fig.  6  where  the  second  harmonic  is  only  20  dB  below  the 
fundamental. 
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Figure  5  -  Power  mode  image  of  a  canine  kidney  before  (left)  and  after  (right)  IV  injection  of 
MfoC- 115  contrast  agent. 
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Figure  6  -  Second  harmonic  signal  amplitude  relative  to  the  fundamental  as  a  function  of  the 
driving  signal  frequency. 
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These  results  for  both  the  leporine  and  canine  experiment  can  now  be  used  as  guides  for  ex 
vivo  experiments.  The  indication  is  that  shorter  than  anticipated  burst  may  be  effective  in  contrast 
interruption  and  time  with  the  cardiac  cycle  may  indeed  be  important.  Most  importantly,  these 
experiments  have  clearly  indicated  that  contrast  interruption  can  be  performed  transcutaneouslv. 
improving  the  chances  of  success  for  this  project.  This  sets  the  stage  for  the  more  extensive  2nd 
year  evaluations  in  animal  studies. 

Decorrelation  Techniques  and  Combination  Blood  Flow  Measurements 

The  manuscript  concerning  decorrelation  techniques  for  perfusion  estimates,  a  significant 
part  of  the  original  proposal,  has  been  refined  during  the  course  of  this  past  year  and  has  now  been 
published.  Attached  is  a  copy  of  the  publication  for  that  work.  Further  direct  progress  in  this  area 
has  been  delayed  to  some  degree  by  a  transition  in  the  hardware  of  the  Diasonics  Spectra  VST 
which  defeated  our  software  as  originally  implemented  on  the  scanner.  The  research  platform  has 
now  stabilized  and  we  are  in  the  process  of  implementing  that  software  again.  Two  improvements 
in  our  situation  with  regard  to  this  software  have  occurred  recently.  The  first  is  that  the  new 
Diasonics  platform  has  a  removable  hard  disk  drive  which  will  make  the  storage  and  transfer  of 
data  more  convenient  and  also  ease  the  transfer  of  new  software  such  as  our  decorrelation  code. 
Secondly,  during  this  downtime  we  have  made  what  may  be  a  significant  use  of  the  decorrelation 
idea  in  another  area  which  impacts  the  technique  which  combines  Doppler  power  and  velocity  data. 
As  was  indicated  in  our  original  proposal,  this  work  is  being  performed  in  conjunction  with  that  of 
Dr.  Paul  Carson  on  3D  breast  imaging.  Even  though  3D  imaging  is  not  a  specific  part  of  this 
grant,  the  idea  was  to  utilize  the  same  data  acquisition  system  so  comparison  could  be  made 
between  techniques.  In  the  course  of  conversations  concerning  the  difficulties  associated  with 
breast  scanning,  it  was  suggested  that  the  elimination  of  encoding  hardware  would  be  desirable  and 
helpful  in  terms  of  examining  images  acquired  from  a  variety  of  scanners.  It  was  suggested  that 
one  might  be  able  to  use  the  decorrelation  of  signal  in  the  elevational  direction  to  determine  the 
separation  between  slices  which  were  obtained  using  a  manual  scan.  The  following  is  a 
description  of  the  techniques  which  is  taken  from  a  preliminary  patent  application. 

The  approach  being  taken  by  most  ultrasound  companies  is  to  build  a  (usually  hand-held) 
motorized  scan  head  specifically  for  3-D  imaging.  This  approach  has  a  disadvantage  of  high  cost, 
low  reliability  and  limitation  on  the  number  of  scan  heads  available  for  3-D.  Often  coupling  paths 
are  employed  with  their  inherent  attenuation  and  reverberation  problems.  Also,  the  scan  rate  for 
optimal  color  flow  imaging  is  a  complex  function  of  all  the  Doppler  time  constants,  frame  rate  and 
the  particular  vascular  signal  levels.  In  regular  imaging  with  a  given  scan  head,  the  operator  or 
observing  radiologist  might  want  to  capture  in  3-D  images  observed  with  a  particular  scanner  using 
an  arcing  or  linear  sweep.  If  relative  3-D  image  position  can  be  obtained  %  processing  the  series 
of  images  or  data  used  to  produce  the  images,  a  great  new  flexibility  and  efficiency  will  be  achieve. 
Major  3-D  capabilities  will  only  be  required  on  a  central  workstation  which  can  record  a  single 
series  of  images  and  display  them  in  their  proper  position  in  a  stacked  slice  3-D  display.  The 
radiologist,  technologist  at  the  scanner,  or  other  trained  observer  will  view  the  slices  sequentially  in 
their  proper  position  on  the  display.  They  will  select  the  most  appropriate  2,  3  or  4-D  display 
method  and  views  to  demonstrate  the  noted  information  in  the  most  informative  way  for  surgeons 
and  other  referring  physicians  or  to  obtain  the  desired  quantitative  measurements. 

Correlation  techniques  have  been  used  extensively  to  estimate  tissue  motion  in  applications 
such  as  blood  flow  (Adler,  Rubin  et  al.  1995,  Appendix  1)  and  imaging  tissue  elastic  properties 
(Adler,  Rubin  et  al.  1990;  Ophir,  Cespedes  et  al.  1991;  Chen,  Jenkins  et  al.  1992;  ODonnell, 
Skovoroda  et  al.  1994).  Some  commercial  ultrasound  systems  now  use  correlative  algorithms  for 
their  color  flow  imaging  rather  than  various  forms  of  Doppler  processing.  We  are  suggesting  the 
use  of  correlation  to  perform  the  opposite  task  namely  assume  the  tissue  to  be  stationary  and  track 
the  movement  of  the  transducer.  Trahey  et  al  (1986)  used  correlation  of  speckle  to  determine  the 
amount  of  lateral  translation  of  a  phased  array  required  to  achieve  statistically  independent 
information  for  speckle  reduction  by  image  compounding.  Their  experimental  results  show 
correlation  curves  as  a  function  of  translation  distance  which  are  smoothly  varying,  suggesting  that 
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measiiring  the  rate  of  decorrelation  and  knowing  the  point  spread  function  for  the  aperture,  one  can 
use  the  same  speckle  statistics  to  estimate  the  actual  transducer  motion.  The  correlation  techniques 
which  we  have  developed  provide  the  mechanism  for  slice  positioning  and  involve  the  adaptation 
of  some  software  already  developed. 

In  order  to  test  the  concept  of  image-based  slice  positioning,  images  of  a  rather 
inhomogeneous  contrast  detail  phantom  were  made.  In  all  cases  the  images  were  post-processing 
in  a  modular  code  under  the  AVS  software  package  of  Advanced  Visualization  Systems.  A  full 
description  of  the  mathematical  process  is  given  in  the  attached  manuscripts  of  Chen  et  al.  Images 
were  collected  from  an  ultrasound  scanner  and  read  into  the  workstation  memory.  The  RGB 
images  are  then  processed  to  select  one  channel  for  B-mode  and  then  the  3-D  data  set  is  sliced  and 
a  single  2-D  plane  displayed  to  select  regions  of  interest  (ROIs)  to  be  processed  for  determining  the 
slice  separation.  The  pixels  contained  in  the  ROI  in  each  image  are  converted  to  a  ID  vector  and  the 
process  repeated  for  each  image.  Each  of  these  vectors  are  then  combined  in  a  2-D  vector.  These 
data  are  then  subsampled  to  determine  the  how  many  of  the  image  ROIs  will  be  used  to  compute 
the  position  of  each  image.  The  group  of  image  ROIs  is  then  andyzed  to  compute  the  correlation 
between  successive  ROIs,  i.e.  ROI#l  correlation  to  ROI#2,  ROI#2  to  ROMS,  etc.  for  a  one  step 
correlation  value.  The  process  is  repeated  for  two,  three,  etc.  step  correlations.  For  the  hand- 
scanned  images,  ROIs  from  groups  of  ten  images  were  used  to  determine  the  slice  separation  for 
the  center  two  slices,  e.g.,  to  position  image  6  with  respect  to  5  use  information  from  images  0-9, 
for  image  7  after  6  use  image  1-10,  etc.  This  does  assume  a  piece-wise  smooth  motion  for  the 
scanhead,  which  should  be  achievable  and  is  necessary  for  any  good  continuous  power  mode 
imaging.  In  each  case,  the  correlation  curve  was  fit  using  a  least  squares  approach  to  the  Gaussian 
functiongiving  the  relationship  between  step  size,  point  spread  function  and  correlation  value.  This 
information  is  then  used  to  position  the  images  appropriately  in  space  for  further  display. 

Two  sets  of  images  were  obtained  in  the  initial  trial.  A  measured  correlation  curve  was  first 
determined  for  mechanical  stepping  with  the  breast  scanner  using  a  3.5  MHz  curved  linear  array. 
Using  several  regions  of  interest  at  the  same  depth  range  and  a  least  squares  method  for  a  Gaussian 

curve-fitting,  the  step  size  was  estimated  for  several  regions  of  interest  to  be  5y  =0.17  ±  0.02  mm, 
compared  to  the  actual  step  size  used  of  0.175  mm.  In  this  case  the  decorrelation  curve  was 
measured  across  a  number  of  images  where  it  could  be  assumed  that  the  step  size  was  constant  for 
the  mechanical  stepping.  A  hand  scan  was  then  performed  where  images  at  32  fps  were  stored  in  a 
120  frame  cine  loop.  Results  given  in  Fig.  7  show  the  appearance  of  the  echogenic  cone  is  not 
very  different  for  this  short  scan  when  using  either  an  equal  slice  separation  assuming  uniform 
hand  scanning  rate  or  that  calculated  using  the  image  decorrelation  rate.  However,  the  lower 
graphs  illustrate  the  decorrelation  does  track  smoothly  the  nonuniform  motion  of  the  slow,  hand 
scanning  currently  used  for  3-D  power  mode  imaging.  Black  stripes  in  the  subtraction  image 
correspond  to  zeros  in  the  integral  displacement  error  of  the  uniform  stepping  as  plotted  in  the 
graph  on  the  right. 

These  techniques  can  be  used  to  correct  positioning  errors  for  hand-scanned  image 
acquisitions  in  the  breast.  This  could  have  an  important  impact  in  the  research  being  performed 
here  since  regional  perfusion  measurements  will  rely  on  know  the  position  of  the  scanhead  when 
combining  Doppler  power  and  velocity  information  together  as  an  estimate  of  perfusion.  The 
image-based  slice  positioning  techniques  will  allow  more  examinations  to  be  performed  in  some  of 
the  more  interesting  cases,  i.e.  the  smaller  dense  breast  which  is  a  particular  difficulty  in 
mammography. 
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Figure  7  -  Left  to  right  top  images;  1)  C-scan,  longitudinal  view  of  a  cone  object  from  a  hand- 
scanned,  3.5  MHz,  curved  linear  array  assuming  equal  step  size;  2)  The  same  view  with  actual  step 
size  as  estimated  by  the  speckle  decorrelation;  3)  The  subtraction  between  these  two  views 
(expanded  horizontally).  For  comparison,  image  1  was  scaled  to  the  same  lateral  size  as  image  2. 
The  left  graph  below  shows  the  step  size  change  with  step  number  as  calculated  from  the 
decorrelation,  while  the  right  graph  shows  the  integrated  displacement  change  with  step  number. 

CONCLUSIONS 

The  results  of  the  research  thus  far  indicates  that  contrast  interruption  is  possible  in  vivo 
using  acoustic  fields  and  such  interruption  results  in  a  negative  contrast  bolus  which  is  well 
resolved  in  time.  Accomplishing  the  transcutaneous  interruption  of  contrast  further  defines  the 
circumstances  which  are  relevant  to  use  in  a  more  complete  study  of  the  process  in  vitro.  These 
studies  should  allow  us  to  reduce  the  acoustic  amplitude  required  to  generate  the  bolus  and  to 
control  the  time  course  of  the  interruption. 

The  development  of  image-based  slice  positioning  should  provide  a  more  accurate 
assessment  of  the  tissue  perfusion  using  the  estimation  techniques  proposed  in  this  research. 
Perhaps  equally  important  is  the  potential  for  these  techniques  to  have  an  impact  on  the  overall  use 
of  3D  ultrasound  to  map  the  breast  vasculature.  Therefore  although  this  was  not  a  specific  goal  of 
this  research,  the  significance  of  this  work  may  reach  beyond  the  concept  originally  proposed  in 
terms  of  decorrelation  measures  and  the  combination  of  Doppler  power  and  velocity  information  to 
estimate  perfusion. 
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Abstract — Imaging  of  blood  flow  perfusion  is  an  area  of  significant  medical  interest.  Recently,  the  advantages 
of  using  the  total  integrated  Doppler  power  spectrum  as  the  parameter  that  is  encoded  in  color  has  been 
shown  to  result  in  an  approximately  threefold  increase  in  flow  sensitivity,  a  relative  insensitivity  to  acquisi¬ 
tion  angle  and  lack  of  aliasing.  We  have  taken  this  mode  a  step  further  and  demonstrated  the  potential  for 
quantifying  blood  flow  using  correlation-based  algorithms  applied  to  the  power  signal.  We  show  that  ^(r) 
=  ^(0)e  T  >  0,  where  ^(r)  is  the  two-time  correlation  of  the  fluctuation  in  the  power  signal,  and  v  is 
the  specific  flow  (reciprocal  of  mean  transit  time).  Scans  of  a  dog’s  blood,  pumped  at  a  constant  rate 
through  gum  rubber  tubing,  were  obtained  using  a  Diasonics  Spectra  10-MHz  linear  array  transducer  at 
standard  range-gated  spectral  mode  (PRF  =  1400  Hz,  wall  filter  =  50  Hz,  sample  gate  =  1.5  mm).  A  fixed 
Doppler  angle  of  68°  was  used.  Five  different  flow  rates  were  tested,  and  the  velocities  determined  by  power 
decorrelation  were  compared  to  the  mean  velocities  calculated  from  the  Doppler  shifts  by  linear  regression 
{R^  =  0.987).  We  believe  the  results  are  very  encouraging  for  using  power  decorrelation  in  perfusion 
evaluation. 

Key  Words:  Decorrelation,  Power  Doppler,  Ultrasound, 


INTRODUCTION 

A  variety  of  physiological  processes,  reflecting  both 
normal  and  pathologic  states,  can  be  characterized  by 
a  relative  increase  in  local  soft  tissue  perfusion,  soft 
tissue  inflammation  being  an  example.  It  is  not  surpris¬ 
ing,  therefore,  that  there  has  been  considerable  interest 
and  expense  devoted  to  producing  quantitative  maps 
of  such  information.  Technological  improvements  in 
ultrasound  imaging  during  the  past  5  years  or  so  have 
shown  this  modality  to  be  capable  of  depicting  blood 
flow  without  the  need  to  administer  intravascular  con¬ 
trast  agents,  although  flow  sensitivity  is  clearly  im¬ 
proved  when  such  agents  are  utilized.  Doppler  ultra¬ 
sound,  using  mean  frequency  maps,  is  a  well-devel¬ 
oped  and  commercially  available  technique  which  can 
noninvasively  depict  flow  in  large  vessels,  provided 
adequate  acoustic  access  exists.  Alternatively,  the  rela¬ 
tively  poor  signal-to-noise  characteristics  of  mean-fre¬ 
quency  maps  has  limited  this  application  as  a  means 


to  fully  and  inexpensively  assess  many  abnormal  flow 
states  commonly  encountered  in  medical  practice. 

The  development  of  a  more  sensitive  ultrasonic 
measurement  of  blood  flow  has  been  proposed  by 
Dymling  et  al.  (1991)  in  which  the  first  frequency 
moment  of  the  backscattered  power  spectrum  has  been 
shown  to  provide  an  estimate  of  flow  in  a  randomly 
oriented  microvascular  bed.  This  latter  property  allows 
one  to  derive  an  explicit  relationship  between  the  first 
moment  and  the  local  tissue  perfusion.  It  is  not  yet 
clear,  however,  as  to  the  general  validity  of  these  as¬ 
sumptions  regarding  the  morphologic  nature  of  the 
capillaries  in  most  tissues.  Alternatively,  the  integrated 
power  spectrum  is  available  on  many  ultrasonic  im¬ 
agers  and,  as  these  authors  point  out,  this  measurement 
is  directly  proportional  to  the  total  number  of  scatterers 
in  the  insonified  volume  (Dymling  et  al.  1991).  The 
nature  of  this  proportionality  may  be  quite  complex, 
however,  being  influenced  by  such  things  as  attenua¬ 
tion,  beam  pattern,  Doppler  wall  filter,  clutter  cancella¬ 
tion  and  frequency  (Eriksson  et  al.  1991 ;  Mo  and  Cob- 
bold  1992). 

Rubin  et  al.  ( 1994)  have  demonstrated  the  utility 
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of  color-encoded  maps  of  integrated  power  on  two 
commercially  available  systems,  referring  to  it  as 
power  Doppler  imaging  and,  more  importantly,  they 
have  demonstrated  a  significant  improvement  in  sig- 
nal-to-noise,  relative  to  mean-frequency  maps,  in  de¬ 
picting  flow.  One  example  consists  of  displaying  the 
soft-tissue  blush  in  the  renal  cortex,  which  typically 
requires  administration  of  an  intravascular  contrast 
agent.  This  particular  feature  is  currently  being  made 
commercially  available  by  a  number  of  ultrasound  ven¬ 
dors  as  a  means  of  increasing  flow  sensitivity  and  more 
accurately  representing  vessel  morphology. 

Inherent  limitations  of  power  Doppler,  besides  its 
significant  motion  sensitivity,  relate  to  more  funda¬ 
mental  issues.  An  estimate  of  the  total  number  of  mov¬ 
ing  scatterers  is  being  displayed  rather  than  the  more 
physiologically  relevant  parameter,  the  blood  flow  it¬ 
self.  The  total  power  displayed  is  not  entirely  velocity 
independent  due  to  the  effects  of  the  soft-tissue  clutter 
canceler,  which  modulates  the  amplitude  of  low  veloc¬ 
ity  components,  thereby  introducing  a  frequency  de¬ 
pendence.  Despite  this  limitation,  we  show,  in  the  cur¬ 
rent  work,  that  the  inherent  stochastic  nature  of  such 
estimates,  prior  to  any  temporal  filtering  {i.e.,  persis¬ 
tence),  can  provide  additional  information.  Such  esti¬ 
mates  have  been  made  in  the  past  based  on  pure  statisti¬ 
cal  fluctuations  in  either  pressure  amplitude  in  A  lines 
or  envelope  variations  from  backscattered  ultrasound 
(Atkinson  and  Berry  1974)  and  using  scattering  theory 
arguments  applied  to  a  Doppler  signal  (Mo  and  Cob- 
bold  1992).  In  our  case,  we  will  estimate  local  tissue 
perfusion  from  the  nontemporally  averaged  data  by 
solving  a  stochastic  differential  equation  derived  from 
a  flow  model  with  noise.  Similar  relationships  are  well- 
known  from  linear  response  theory  and  relate  the  fluc¬ 
tuations  that  drive  a  system  to  the  systematic  response 
of  the  system  {i.e.,  a  transport  coefficient  in  the  case 
of  hydrodynamic  phenomena)  via  the  fluctuation  dissi¬ 
pation  theorem  (Mori  and  Fujisaka  1973;  Papoulis 
1965;  Zwanzig  1960).  Inasmuch  as  such  a  mathemati¬ 
cal  approach  to  local  biologic  phenomena  is  not  partic¬ 
ularly  well-known,  the  relevant  theoretical  considera¬ 
tions  are  developed  in  some  detail  below.  It  is  shown 
that  the  rate  of  decorrelation  within  a  region-of-interest 
directly  relates  to  specific  flow,  as  suggested  below.  To 
illustrate  the  technique,  a  simple  flow  phantom  under 
conditions  of  constant  input  is  studied.  Mean  velocity 
(proportional  to  flow  in  this  geometry)  is  directly  com¬ 
pared  to  estimated  velocity  based  on  known  sample 
volume  and  estimated  transit  time  obtained  from  decor¬ 
relation. 

THEORY 

The  principal  assumption  in  power  Doppler  im¬ 
aging  pertains  to  the  direct  relationship  between  total 
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power  within  a  voxel  (or  region  of  interest)  and  the 
mean  number  of  scatterers  within  that  voxel.  This  rela¬ 
tionship  holds  provided  that  multiple  scattering  events 
are  negligible  for  the  range  of  fractional  blood  volumes 
considered,  and  the  flow  is  not  turbulent.  Typically 
power  Doppler  images  are  displayed  with  large 
amounts  of  temporal  averaging  with  different  frame 
weightings  to  diminish  the  effects  of  soft  tissue  motion 
(Rubin  et  al.  1994).  Such  maneuvers  diminish  the  sta¬ 
tistical  fluctuations  in  the  power  display.  To  restore 
these  fluctuations,  the  temporal  averaging  is  removed. 
Thus,  for  each  frame  within  the  nontemporally  aver¬ 
aged  “power”  image,  a  number  N  is  derived,  repre¬ 
senting  the  instantaneous  estimate  of  the  number  of 
scatterers  entering  the  tissue  volume,  V,  per  unit  time. 
This  number,  in  general,  is  subject  to  some  variation 
from  frame-to-frame  due  to  inherent  noise  in  the  rate 
of  scattering,  as  well  as  the  acquisition  and  signal  pro¬ 
cessing  system  itself.  These  factors  include  attenuation 
and  beam  geometry  as  mentioned  earlier.  However, 
the  final  specific  flow  is  obtained  from  a  normalized 
expression,  canceling  many  of  these  effects.  A  stochas¬ 
tic  equation  is  derived  which  contains  information  re¬ 
lated  to  the  rate  at  which  scatterers  enter  the  voxel, 
provided  that  certain  specific  assumptions  are  valid 
(see  below).  For  simplicity,  we  consider  a  constant 
input. 

We  suppose  that  the  input  flow  is  not  random  and 
depends  only  on  the  number  of  scatterers  (RBCs  in 
this  case)  and  flow  per  unit  mass,  F.  In  general,  since 
the  instantaneous  estimate  of  N  is  subject  to  noise, 
the  actual  and  estimated  inputs  are  not  equivalent.  To 
describe  the  time  dependence  of  this  random  variable, 
therefore,  it  is  necessary  to  write  an  appropriate  sto¬ 
chastic  differential  equation.  In  this  case  we  assume 
the  simple  form: 

dN/dt  =  F{Ca  -  mN/V)  +/(f)  (1) 

where  m  =  mass  per  scatterer,  F  =  flow  per  unit  mass, 
V  —  volume  of  interest,  N  =  instantaneous  number  of 
scatterers  in  V,f(t)  =  stochastic  driving  force,  Ca  — 
arterial  mass  concentration  of  scatterers. 

The  time  rate  of  change  in  N  depends  on  the  sys¬ 
tematic  term  (Le.,  due  to  flow)  and  a  noise  term,/(r), 
characterizing  the  stochastic  nature  of  the  system.  This 
so-called  Langevin  equation  is  well-known  within  the 
statistical  mechanics  literature  (Mori  and  Fujisaka 
1973;  Papoulis  1965).  Several  formal  properties  of 
this  relation  may  be  derived;  these  may  be  regarded 
as  a  subset  of  the  more  general  Langevin  approach 
(Mori  and  Fujisaka  1973).  Upon  averaging  eqn  (1) 
over  all  realizations  of  A,  recognizing  that  under  steady 
flow  conditions  (N)  is  constant,  results  in 
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d{N)ldt  -  0.  (2) 

Farther,  by  conservation  of  mass: 

FC^  =  Fm{N)fV 
or 

(N)  =  QWm,  (3) 

The  mean  number  of  scatterers  contributing  to  the 

power  image  reflects  mean  transport  of  RBCs  into  V. 
From  eqs  ( 1 )  and  (3),  it  follows  that: 

{/(O)  =  0.  (4) 


6N(t)  =  e-'^WiO)  +  f  dTe-''^-'^f{T)  (10) 

do 

which,  after  forming  #(r)  =  {6N{t  +  t )6N{t)),  re¬ 
sults  in 

$(r)  =  [$(0)  -  +  #(0)^-''''-'  '  (11) 

in  which  \  \  denotes  absolute  value.  For  stationarity 

to  be  satisfied,  we  require  that 

a  =  2i/$(0)  (12) 

so  that  eqn  (11)  simplifies  to 

#(t)  =  *(0)c-"L  r>0  (13) 


The  fluctuation  in  N 


6N  =  N  -  (N)  (5) 

satisfies  a  somewhat  simplified  expression 


d6N 

dt 


-'y6N+fit)=-u6N+f(,t)  (6) 


where 


1/  =  mFIV  —  specific  flow.  (7) 

The  remaining  points  to  consider  deal  specifically  with 
the  statistical  properties  of  this  system.  For  this  pur¬ 
pose,  a  new  function  is  introduced 

^(r)  =  {6Nit  +  T)6N{t))  (8) 

which  is  the  two-time  correlation  function  for  the  ex¬ 
pected  number  of  scatterers  in  V.  Three  additional 
properties  must  first  be  developed.  Stationarity  is  al¬ 
ready  inherent  in  eqn  (8)  and  indicates  the  indepen¬ 
dence  of  ^(r)  on  the  initial  time  t.  The  second  prop¬ 
erty,  termed  causality,  indicates  that/(r)  and  6N(t*) 
are  not  correlated  when  t'  <  t.  The  final  statement 
refers  to  the  basic  Markov  nature  of  eqn  (6),  which 
when  combined  with  the  assumption  of  stationarity, 
results  in  the  so-called  fluctuation  dissipation  theorem: 

{mf{F))  ^  a6(t  -  F)  (9) 

in  which  it  is  apparent  that/(t)  is  a  white  noise  process 
and  cr  is  a  constant.  In  terms  of  these  assumptions,  a 
simple  expression  for  $(t)  will  be  derived.  Solution 
of  eqn  (6)  results  in 


which  is  to  say  that  the  decay  of  the  correlation  func¬ 
tion  itself  gives  an  estimate  of  the  specific  flow. 

EXPERIMENTAL  METHODS 

Experiment  1 

Preliminary  experiments  were  performed  to  test 
the  accuracy  of  the  decorrelation  method  in  the  mea¬ 
surement  of  perfusion.  In  a  simple  flow  experiment, 
gum  rubber  tubing  with  an  internal  diameter  1  cm  was 
positioned  so  that  the  angle  of  incidence  was  0°  be¬ 
tween  the  flow  stream  and  imaging  field  of  a  Diasonics 
Spectra  VST  using  a  10-MHz  probe.  This  geometry 
maximized  the  Doppler  frequency  shift  for  any  flow 
velocity.  The  flowing  media  was  a  mixture  of  com 
starch  (2  g/L)  in  water  circulated  by  a  peristaltic  pump 
with  the  fluid  reservoir  acting  as  a  capacitance  chamber 
to  produce  a  constant  flow  of  10  to  20  mL/min  through 
the  tube.  Knowing  the  size  of  the  pixels  used  in  dis¬ 
playing  the  power  mode  color  and  the  input  volume 
flow  rate  to  the  tube,  it  was  estimated,  assuming  a 
parabolic  velocity  profile  across  the  tube,  how  rapidly 
each  pixel  would  have  to  be  sampled  to  observe  the 
decorrelation  time  at  various  positions  within  the  tube. 
The  rate  at  which  the  power  color  display  is  updated 
at  any  given  pixel  corresponded  to  the  frame  rate  of 
the  ultrasound  imager.  At  this  rate,  all  but  the  slowest 
flows  essentially  decorrelated  instantly.  To  lengthen 
the  decorrelation  time,  effective  pixel  size  was  in¬ 
creased  by  summing  pixels  along  the  direction  of  flow 
and  using  the  total  power  in  these  pixels  as  if  they 
represented  one  sampling  site.  Hence,  the  total  power 
in  n  contiguous  pixels  along  the  flow  direction  effec¬ 
tively  lengthens  the  region-of-interest  and  permits 
measurement  of  decorrelation  rates  that  are  n  times 
faster  than  before  {Le,,  the  transit  time  is  increased  by 
a  factor  of  n). 


496 


Ultrasound  in  Medicine  and  Biology 

Using  software  made  available  by  Diasonics  for 
experimental  use  on  their  Spectra  machine,  a  buffer  is 
filled  with  sequential  frames  of  nontemporally  aver¬ 
aged  power  mode  data.  This  data  is  then  descrambled 
and  loaded  into  an  image  processing  software  package 
( Application  Visualization  System  [  AVS  ] )  in  the  form 
of  a  three-dimensional  (3D)  (2  spatial,  1  temporal) 
scalar  byte  field.  Each  byte  represents  a  decibel-scale 
value  of  the  Doppler  signal  amplitude  at  that  location 
For  each  set  of  data,  a  region  of  interest  is  selected 
within  an  area  of  constant  flow,  i.e.,  a  specified  distance 
from  the  center  of  the  tube.  The  pixels  within  this 
region  are  written  to  one  line  of  a  two-dimensional 
(2D)  field.  Each  line  of  this  field  is  sequentially  filled 
from  the  same  ROI  (region-of-interest)  of  each  frame 
in  the  original  3D  field. 

Each  point  in  this  field  is  then  converted  to  a 
linear  scale  and  the  global  mean  is  subtracted. 


1 


M-l 


Fix,  t)  =  -  M  ^ 

m  —0 


Hx.nt )/ 106.67 


where  P(x,  r)  is  the  linear  zero  mean  power  value, 
I{x  f)  is  the  raw  decibel  data,  M  is  the  number  o 
frames  and  L  is  the  size  of  the  ROI  minus  ;  -  1, 
where  7  is  the  number  of  pixels  summed  along  the  flow 
direction  to  “lengthen”  a  pixel.  The  factor  106.67  is 
used  to  convert  the  0  to  255  log  scale  for  intensity  to 
a  linear  scale. 

The  correlation  is  then  calculated  using  the  ex¬ 
pression 


corr(m)  =  -^  S  IL  PU,  n)  X  PU,  n  +  m) 

M  -  m  /=o 

where  corr  (m)  is  the  time  correlation  between  ROIs, 
and  m  is  the  frame  separation. 

The  normalized  correlation  is; 


Corr  (m) 


corr  (m) 
corr  (0) 


Experiment  2  ,  i 

In  a  second  experiment,  the  spectral  Doppler 

mode  was  used  to  obtain  data  at  a  much  faster  rate 
than  was  possible  in  the  power  mode.  This  allowed 
measurement  of  faster  flows  where  the  Doppler  shift 
was  more  easily  detected  in  addition  to  the  greater  flow 
sensitivity  of  the  spectra  Doppler  mode  itself.  In  this 
case,  fresh  canine  blood  was  used  as  the  flowing  me¬ 
dium  through  a  length  of  0.6-cm-diameter  dialysis  tub¬ 
ing.  The  pumping  system  was  unchanged  and  the 
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Doppler  spectral  data  was  acquired  from  the  center  of 
the  tube  using  a  fixed  angle  of  68“  at  various  flow 
velocities  as  measured  from  the  Doppler  spectrum,  lo 
analyze  the  spectral  data,  the  digital  data  was  again 
dumped  from  the  Diasonics  Spectra  and  analyzed  using 
AVS  The  analysis  consisted  of  converting  the  digital 
values  coded  in  8  bits  to  the  42-dB  dynamic  range  of 
the  Doppler  spectrum.  In  general,  spectral  Doppler  data 
consists  of  a  series  of  power  spectra  which  are  acquired 
at  each  time  step  and  encoded  in  pixel  intensity.  There¬ 
fore,  at  each  time  point  in  the  data,  the  corresponding 
power  spectrum  was  integrated  to  yield  the  same  result 
for  the  single  Doppler  sample  volume  as  is  encoded 
in  color  at  each  pixel  in  the  power  mode  image.  The 
sampling  interval,  corresponding  to  the  time  windows 
over  which  individual  fast-Fourier  transforms  were 
performed,  is  about  2  ms.  This  is  about  50  times  faster 
than  the  frame  rate  in  the  color  mode,  typically  about 
10  frames/s.  The  integrated  power  was  then  windowed 
and  an  autocorrelation  performed  with  the  mean  sub¬ 
tracted.  The  decorrelation  results  were  then  used  to 
compute  the  rate  at  which  the  blood  was  flowing 
through  the  Doppler  sample  volume,  and  this  result 
was  compared  to  the  mean  flow  velocity  as  measured 
using  the  angle  corrected  Doppler  frequency  shift  ac- 
cording  to 


f  fPif)df 

/  =-p - 

J  Pi  f)df 

where  P(f )  is  the  Doppler  power  spectrum  and  /  is 
the  mean  frequency  shift. 

results 

Experiment  1  ,  •  u  ^  u 

Figure  la  depicts  the  decorrelation  within  the  tube 

flow  at  two  locations.  Assuming  laminar  flow,  which 
is  appropriate  for  these  flow  rates,  the  decorrelation  is 
faster  in  the  center  of  the  tube  than  near  the  walls,  as 
expected.  Using  the  ROI  near  the  tube  wall,  the  number 
of  pixels  treated  as  a  single  sample  region  was  varied 
and  the  effect  of  this  is  presented  in  Fig.  lb.  As  the 
number  of  pixels  is  increased,  effectively  lengthening 
the  ROI  the  decorrelation  time  increases.  This  is  con¬ 
sistent  with  the  notion  that  the  fluid  flowing  through 
this  enlarged  region  would  take  a  long  penod  of  time 
to  pass. 

Experiment  2 

Using  the  spectral  Doppler  data,  the  results  are 
similar  to  those  obtained  by  the  power  mode  imaging. 
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Time  (sec)  Time  (sec) 

(a)  (b) 

Fig.  1.  Results  taken  using  the  power  mode  decorrelation  technique  in  the  tube  flow  system.  In  this  case,  the  tube 
was  positioned  such  that  the  flow  was  directly  toward  the  transducer.  Therefore,  pixels  could  be  averaged  along 
the  flow  direction  (see  Methods).  The  pixel  number  indicates  the  number  of  pixels  used  in  the  average,  (a)  The 
data  indicates  the  effect  on  the  decorrelation  rate  for  the  increased  flow  in  the  center  of  the  vessel  compared  to 
the  flow  near  the  walls,  (b)  A  comparison  of  the  effect  of  the  number  of  pixels  used  in  the  average. 


An  example  of  the  integrated  power  spectra  versus 
time  is  given  in  Fig.  2a.  Note  the  quasiperiodicity  in 
the  integrated  power  related  to  the  speckle  pattern  of 
the  finite  aperture  used  in  the  ultrasound  imaging.  Fig¬ 
ure  2b  is  a  plot  of  the  decorrelation-calculated  flow 
velocity  versus  the  flow  velocity  as  measured  by  the 
Doppler  frequency  shift.  The  conversion  of  the  decor¬ 


relation  time  to  a  velocity  is  based  on  the  transit  of 
the  Doppler  sample  volume  in  the  direction  of  the  flow 
and  the  time  for  decorrelation  to  0.5.  The  slope  of  the 
linear  fit  is  0.91  with  an  =  0.98  indicating  nearly 
a  one-to-one  correspondence  for  the  measured  veloci¬ 
ties.  The  intercept  for  the  fit  is  0.69  mm/s,  which  indi¬ 
cates  the  movement  of  blood  in  and  out  of  the  Doppler 


Time  (msec)  Doppler  Velocity  (mm/sec) 


Fig.  2.  Decorrelation  results  using  spectral  Doppler  data,  (a)  The  integrated  power  spectra  at  each  time  point. 
Note  the  appearance  of  a  periodicity  associated  with  the  movement  of  speckle,  (b)  A  comparison  of  the  velocity 
calculated  by  the  decorrelation  and  the  standard  Doppler  shift  technique.  The  line  (slope  =  0.91,  /?“  =  0.98)  is 

a  linear  fit  of  the  data  points. 
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sample  volume  when  there  is  no  net  flow.  This  situa¬ 
tion  can  occur  in  this  case  due  to  the  large  vessel 
diameter.  It  is  expected  that  flow  in  small  vessels  (arte¬ 
rioles,  venules,  and  capillaries,  if  Doppler  shitt  is  de¬ 
tected)  will  be  more  deterministic  due  to  the  greater 
physical  constraints  on  the  blood  motion  in  the  smaller 
vessels. 

DISCUSSION 

Because  of  the  ubiquitous  relationship  of  perfusion 
abnormalities  to  disease  states  and  their  therapies,  a 
simple  and  accurate  in  vivo  measurement  of  pertusion 
is  highly  desirable.  Techniques  that  measure  perfusion 
employ  computed  tomography  (CT),  magnetic  reso¬ 
nance  imaging  (MRI)  and  nuclear  medicine;  these  typi¬ 
cally  analyze  kinetic  profiles  of  contrast  agents  to  e^mct 
tissue  blood  flow  parameters  (Beaney  et  |  984  Hu¬ 
ang  et  al.  1982;  Kallinowski  et  al.  1989;  ^  ’ 

Miles  et  al.  1993;  Peters  et  al.  1987;  Rosen  et  al.  1990). 
Beside  the  obvious  need  to  inject  contrast  agents,  these 
agents  usually  are  not  confined  to  the  blood  stream^ 
Hence,  difficulties  arise  in  extracting  true  tissue  blood 
flow  from  other  parameters  such  as  leakage  of  the  agents 
into  the  soft  tissue  compartment.  Some  of  this  comp  ex- 
ity  can  be  removed  with  appropriate  kinetic  modeling 
( Adler  1 987 ) ,  but  it  seems  doubtful  that  all  of  the  multi¬ 
ple  variables  involved  can  be  accounted  for. 

Ultrasound  has  some  very  attractive  features  or 
perfusion  measurements.  Since  red  blood  cells,  an  in¬ 
herently  intravascular  nondiffusible  substance,  are  the 
primary  source  of  ultrasound  backscatter  in  blood,  only 
perfusion  is  effectively  being  measured.  Second,  since 
most  ultrasound  contrast  agents  are  microbubble  based, 
they  not  only  improve  signal-to-noise  ratio  but  are  non¬ 
diffusible  as  well,  unlike  typical  x-ray  or  magnetic 
resonance  contrast  agents.  These  ultrasound  agents 
have  already  provided  clinically  useful  information  tor 
perfusion  of  such  tissues  as  the 

strong  1986;  Kaul  1990;  Wilson  et  al.  1993).  Third, 
the  nonlinear  oscillation  of  bubbles  makes  it  possi  e 
to  separate  flow  signal  from  background  signals  even 
further  by  imaging  the  second  harmonic  of  the  trans¬ 
mitted  ultrasound  frequency  (Shrope  and  Newhouse 
1993).  The  bubble  response  at  the  second  harmonic  is 
proportionately  much  greater  than  the  relatively  weak 
response  of  static  tissues.  This  tends  to  amplify  flow 
signals  relative  to  soft  tissue  making  pertusion  mea¬ 
surements  much  more  practicable. 

Dymling  et  al.  ( 1991 )  used  the  first  moment  ot 
the  Doppler  power  spectrum  for  calculating  perfusion. 
The  method  employs  CW  Doppler  and  has  shown 
promise,  working  in  a  phantom  and  in  a  volunteer  s 
finger.  The  method  is  somewhat  geometry  dependent, 


and  relies  on  an  isotropic  distribution  of  flows  for  an 
accurate  estimate  of  perfusion;  however,  the  promise 
of  a  contrast-free  perfusion  measure  is  intnguing. 

Power  Doppler  imaging  is  a  new  ultrasonic  tech¬ 
nique  for  displaying  blood  flow  with  significantly  im¬ 
proved  sensitivity  relative  to  conventional  color  Dopp¬ 
ler  imaging.  It  is  essentially  angle  independent,  and  it 
is  not  subject  to  aliasing  artifacts.  Vessel  morphology 
is  well  depicted  using  this  technique,  and  the  microvas- 
cular  bed  in  some  tissues  has  been  shown  to  appear  as 
a  blush.  This  technique  has  shown  considerable  prom¬ 
ise  in  mapping  areas  of  perfusion,  or  its  absence,  in  a 
number  of  clinical  situations  (Rubin  1993,  in  press). 
Detailed  clinical  studies  are  currently  underv^y  to  bet¬ 
ter  define  its  role  (Newman  et  al.  1994).  The  major 
drawbacks  of  the  technique  in  its  current  foim  rel^ate 
to  motion  sensitivity  and  loss  of  velocity-related  infor¬ 
mation.  Further,  in  its  current  state,  a 
to  local  blood  volume  rather  than  the  actual  blood  flow 

is  being  encoded. 

Atkinson  and  Berry  ( 1974)  have  shown  that  the 
fluctuations  of  ultrasound  pressure  amplitude  or  fluc¬ 
tuations  in  the  envelope  derived  from  the  pressure  am¬ 
plitude  can  generate  information  about  motion  through 
the  beam.  Approaching  the  problem  from  a  purely  sta¬ 
tistical-based  scattering  theory  in  contrast  to  the  dy¬ 
namic  approach  we  employed,  they  used  a  time  domain 
cross-correlation  technique  to  derive  “fading  times 
that  are  similar  to  the  decorrelation  times  we  calcu  ate 
from  the  power  in  the  frequency  domain  (Doppler) 
signals.  In  fact,  their  Fig.  lb  looks  remarkably  similar 
to  our  Fig  2a.  This  is  of  further  interest,  since  their 
“fading  time”  is  based  on  motion  transverse  to  the 
sound  beam,  whereas  our  spectra  were  obtained  at  68  , 
thus  containing  both  transverse  and  longitudinal 
speckle  components.  Since  the  properties  of  speckle 
are  different  across  the  beam  as  opposed  to  ^ong  the 
beam  (Burkhardt  1978),  we  are  encouraged  that,  at 
least  in  this  single  intermediate  case,  the  results  were 
somewhat  immune  to  speckle  morphology.  Of  course, 
our  experiment  1,  which  employed  color  power  im¬ 
aging,  only  measured  motion  parallel  to  sound  beam 

*’^'^Alfo  ’ approaching  this 

scattering  theory  arguments.  Mo  and  Cobbold  ( 1992) 
have  suggested  that  the  fluctuations  in  the  amplitude 
of  the  Doppler  signal  relate  to  the  rate  at  which  scatter- 
ers  move  through  a  sampling  region.  This  is  manifested 
as  a  decorrelation  time  analogous  to  ours.  It  is  interest¬ 
ing  that  similar  conclusions  can  be  drawn  regarding 
dynamical  behavior  of  the  power  in  the  Doppler  signal 
whether  starting  from  a  scattering  theory  approach  or 
using  purely  kinetic  arguments  such  as  ours. 

The  techniques  using  fluctuations  in  the  Doppler 
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pow'er  have  some  definite  advantages  over  the  A-mode 
methods.  The  Doppler  methods  can  be  easily  imple- 
rnented,  since  the  Doppler  power  spectrum  is  already 
calculated  during  standard  spectral  analysis,  and  the 
integration  of  these  values  to  obtain  the  total  power 
would  be  trivial.  Second,  the  Doppler  spectrum  is  al¬ 
ready  clutter  canceled,  removing  soft-tissue  motion,  a 
significant  source  of  noise.  In  fact,  effective  clutter 
cancellation  already  employed  in  obtaining  the  Dopp¬ 
ler  signal  may  be  one  of  the  biggest  advantages  in 
using  the  Doppler  method,  since  the  overwhelming 
signal  from  stationary  tissue  could  confound  any  of 
these  decorrelation  techniques.  While  approached  from 
different  starting  points,  the  two  methods  exhibit  simi¬ 
larities,  and  A-line  decorrelation  techniques  have  al¬ 
ready  been  used  with  success  to  image  blood  flow 
(Bamber  et  al.  1988). 

We  have  demonstrated,  in  two  simple  flow  phan¬ 
toms,  that  the  technique  may  allow  one  to  measure 
local  tissue  perfusion  using  a  correlation-based  esti¬ 
mate  of  mean  transit  time.  In  fact,  the  mean  velocity 
calculated  from  the  normalized  first  moment  of  the 
Doppler  power  spectrum  has  a  one-to-one  correspon¬ 
dence  (slope  =  0.91,  =  0.98,  Fig.  2b)  with  the 

velocity  calculated  from  mean  transit  time  estimated 
from  the  decorrelation  profiles  and  the  known  sample 
volume  size  over  a  broad  range  of  flows  employed  in 
the  second  experiment.  A  critical  feature  in  obtaining 
these  estimates  relates  to  sampling  rate;  undersampling 
can  result  in  a  situation  analogous  to  aliasing  in  con¬ 
ventional  Doppler,  producing  confusing  decorrelation 
spectra.  This  is  not  surprising  given  the  quasiperiodic 
nature  of  the  raw  power  data  (Fig.  2a).  Alternatively, 
at  sufficiently  slow  flow  rates  (expected  in  true  capil¬ 
lary  flow  situations),  simple  exponential  profiles  were 
obtained,  which  resulted  in  unambiguous  estimates  of 
mean  transit  time. 

A  major  assertion  in  this  analysis  is  the  validity 
of  eqn  ( 1 ),  which  requires  that  there  be  no  other  time 
scale  comparable  to  the  mean  transit  time  through  the 
region  of  interest.  Second,  the  nature  of  the  measure¬ 
ments  being  made  are  assumed  to  represent  a  stationary 
process,  which  seems  to  be  a  reasonable  assumption 
(Mo  and  Cobbold  1992).  Finally,  we  have  assumed 
constant  flow  both  locally  as  well  as  temporally.  It  is 
expected  that  the  theory,  as  stated,  can  be  generalized 
to  more  complex  temporal  flow  dependencies.  Further 
it  is  expected  that  in  most  situations  a  region  of  rela¬ 
tively  constant  flow  can  be  selected  so  that  one  may 
avoid  polydispersity  in  the  distribution  of  flows  for  a 
given  region. 

This  point  may  be  moot,  however,  since  the  fluc¬ 
tuations  in  power  that  we  are  detecting  are  almost  cer¬ 
tainly  not  due  to  individual  red  blood  cells  but  are 


caused  by  changes  in  coherent  scattering,  speckle,  as 
groups  of  scatterers  move  through  the  sound  beam. 
The  effect  of  speckle  moving  in  and  out  of  the  Doppler 
sample  volume  has  produced  a  quasiperiodic  spectrum 
even  though  the  flow  itself  is  constant  (Fig.  2a).  This 
fact  alone  may  ultimately  necessitate  a  more  complex 
theoretical  model  beyond  the  simple  exponential  de- 
correlation  demonstrated  herein,  even  for  constant 
flows.  The  entire  process  would  be  further  complicated 
by  the  presence  of  turbulence,  where  multiple  local 
velocities  will  cause  rapid  fluctuations  in  power  due 
changes  in  scatterer  density  or  coherent  interference 
(Angelsen  1980;  Shung  et  al.  1992;  Twersky  1978). 
These  changes  would  not  reflect  actual  flow,  although 
one  could  envision  using  such  rapid,  unexpected  decor¬ 
relations  as  a  measure  of  turbulence.  Such  circum¬ 
stances  would  be  highly  unlikely  in  the  slow  flow, 
small  vessel  situations  that  we  are  interested  in  investi¬ 
gating,  however.  Nevertheless  even  with  these  poten¬ 
tial  problems,  we  have  shown  that  the  two-time  power- 
derived  correlation  within  a  region-of-interest  can  be 
used  to  estimate  mean  flow.  We  currently  plan  to  per¬ 
form  experiments  in  more  physiologic  flow  phantoms 
(e.g.,  excised  kidneys,  etc.)  in  which  quantitative  esti¬ 
mates  of  perfusion  in  randomly  configured  vascular 
beds  can  be  performed. 

We  have  shown  that  the  fluctuations  inherent  in 
deriving  an  estimate  of  mean  scatterer  number,  and  the 
basic  multiframe  nature  of  ultrasonic  imaging,  pro¬ 
vides  a  simple  scheme  to  estimate  local  tissue  perfu¬ 
sion.  The  time  scales  over  which  fluctuations  in  mean 
scatterer  number  decorrelate  may  be  regarded  as  mean 
transit  time  for  passage  of  RBCs  through  the  tissue 
volume  of  interest.  We  have  shown  this  number  to  be 
a  direct  measure  of  specific  flow.  Having  obtained  this 
information,  a  true  perfusion  map  may  be  derived. 
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Abstract 


In  this  paper,  the  correlation  function  of  the  echo  signal  intensities, 
which  is  directly  related  to  the  changes  of  speckle  pattern  on  a  series  of  B- 
mode  images,  is  used  to  determine  the  scan-plane  positions,  motions,  and 
orientations.  It  is  assumed  that  the  statistical  properties  of  the  echo  signals 
follow  those  of  a  complex,  circular  Gaussian  and  the  case  is  considered  of 
diffusely  scattering  tissue  with  many  fine  particles  per  resolution  cell  and  with 
no  phase  distortion.  The  random  processes  involved  in  forming  these  echo 
signals  give  rising  to  an  uncertainty  in  the  estimate  this  correlation  function. 
This  statistical  uncertainty  is  analyzed  using  error  propagation.  Finally,  the 
method  is  applied  to  1-D  array  transducer,  and  initial  results  are  presented 
from  a  tissue-mimicking  phantom.  Experimental  results  are  in  reasonable  good 
agreement  with  the  predictions.  The  method  described  in  this  paper  can  be 
applied  to  a  clinical  set  when  tissue  can  be  modeled  as  a  gelatin  matrix 
containing  many  small  randomly  distributed  point-like  scatterers  with  uniform 
acoustic  properties.  The  current  method  should  serve  as  a  baseline  for  the 
motion  detection  in  general  clinical  applications. 
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INTRODUCTION 

T’u  rvf  thw  naner  is  to  studv  the  scan-planc  motions  through  tissue  and  the 

The  purpose  of  this  paper  nf  sneckle  correlation  in  ultrasonic  images, 

estimate  of  their  relative  positions  y  .  J  the  correlation  of  speckle 

Conventionally,  the  second-order  -  -  g  tracking  the  relative  motion  of  blood 

patterns  at  been  sfudied  for  the  purpose  of  the 

or  other  moving  tissue  [1].  Speck  ,  .  spatial  and  frequency  compounding  [2, 

speckle  reductton  ta^  —  SeThe  averaging  of  images  in  which  the 

3,  4,  and  5].  These  specKie  speckle  correlation  between 

mg?s“o?S»  Se“Sve  detection  of 

resolution  cell  of  the  transaucer.  wucu  ^  chanee  producing  a  change  in  the 

domain  echo  signals  at  one  distance  alo  g  iniases  We  will  analyze  the  second- 

spcckie  patterns  observed  at  that  location  on  function  of  time  or 

“pS  at  two  posiUons  in  space,  number  of  frame  or 

In  this  paper,  it  is  assumed  that  the  ^latistical  pro^ri^of  echo 
complex,  circular  Gaussian  and  *e  c^c^s  con  The  method,  in  which  the 

fine  particles  per  resolution  cell  PI  and  ph^^^  deteimine  the  relative  positions  of  scan 

^;aTe‘frdto°rfiotrifprpposed.  Then,.he^s«^^^^^^ 

!i“rr"4‘ »n“and  thelniW  Ste  rpLnted  from  testing  ^  tissne-tnimicking  phantom. 

I.  THEORY 

raftSn1)1ase?:SSl?rdt"mr 

comltton  of  the  echo  signal  intensities  between  two  specific  positions  separated  by  dr ,  e.g., 
r  and  r  -I-  4? ,  is  defined  as 


R{Ar)  =-<  Iir)I{r  +  Ar)> 


(1) 


where  7(r)  =  |r/(rf  and  /(r  +  df )  .  |l/(f  k  drf  are  the  intensities  of  the  echo  signals  at 
specific  two  points  r  and  r+df.  «?)  and  I/(r  +  dr)  are  the  complex  echo  signals.  <...> 

normal,  zero  mean,  Gaussian  random 
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variables,  and  assuming  that  the  real  and  imaginary  parts  of  the  complex  echo  signal  f/(r), 
[/^(r)  and  C/,  (r),  are  uncorrelated,  e.g. 


{U,  {r)Ui{r  +  Ar))  =  {U^ir  +  Ar)Ui  (r))  =  0, 

we  have  [7] 

R{Ar)  =  2(U ? {r)^(Urir  +  Ar)^  +  4{t/ ^ (r)[7^(r  +  Ar))^  +  2(uf  (.r)^(Ui  (r  +  Ar)j . 
Using  the  symmetry  of  the  correlation  function  of  the  complex  echo  signals  f/(r). 


Finally,  we  have 
or 


p{Ar)  = 


(u(r)U*  {r  +  Ar)) 


</> 


R{Ar)  =<  /  >^  (1  +  11p( Ar )lf ) 


iiP(AF)r=^-i 


(2) 

(3) 

(3') 


IS 


where  <  /  >  is  the  average  echo  signal  intensity  over  the  region  of  interest.  Here  llp(Ar 
called  the  normalized  correlation  of  the  echo  signal  intensities  and  is  denoted  by  r(  Ar) . 

Equations  (3)  and  (3')  show  that  there  is  a  simple  relationship  between  the  correlation 
function  of  the  complex  echo  signals  and  the  correlation  function  of  the  signal  intensities,  which 
is  directly  related  to  the  speckle  on  B-mode  images. 


b.  Applying  for  clinical  scanners 

The  ideal  instruments  for  the  3-D  ultrasonic  imaging  currently  are  array  scanners,  here  a 
1-D  linear  array  will  be  considered.  To  further  simplify  our  analysis,  we  confine  our 
measurements  in  the  focal  zone  or  Fraunhofer  region  of  the  away  transducer.  Applying  the 
Huygens-Freanel  principle,  the  radiating  surface  of  the  transducer  as  a  collection  of  point 
sources  that  oscillate  sinusoidally  is  considered.  If  the  sources  he  entirely  within  a  simple  plane 
(x',y')  as  shown  in  Figure  1,  the  radiating  surface  is  defined  by  an  aperture  distribution 

function  [8] 


a(x',y')' 
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where  ix',y')  is  a  point  on  the  transducer  surface.  W  is  the  width  of  array  elements,  d  is  the 
center-to-center  distance  between  array  elements,  D  is  the  array  transducer  dimension,  h  is  the 
height  of  array  elements  in  the  elevation  direction,  and  functions 
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rect{x)  =  ^ 
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When  the  pulsed  transducer  is  used  to  insonify  a  medium  containing  randomly 
distributed  fine  particles,  the  complex  echo  signal  is  given  by  [8] 


likz 

U{r)~C-^\^{x,y) 


(5) 


where  k^lKlk  is  the  ultrasound  wave  number,  is  ultrasound  wave  length  at  its  central 
treauencv  z  is  a  distance  along  beam  axis  from  the  surface  of  the  transducer  to  the  field  porn  , 
and  h,  J’y)  is  called  the  directivity  function  and  is  specified  entirely  in  the 
Se  regton  of  the  interest )  when  the  assuiytion  of  local  plane-wave  exists.  For  a  1-D  array 
transducer,  the  directivity  function  is  given  by  [8J 


+00 


^sinc( 


f=-o 


wx  ,  .  ,  Dx  iD. 

— )sinc(- - j) 


sincA 

AoZ 


(6) 


H.rp  the  motion  of  scan-plane  in  the  perpendicular  direction  related  to  scan  plane  (called 
given  by 

r(Ay,.)  =  (sincVo  ®  sinc^J^ 

where  0  stands  for  a  convolution  operator.  Ay,  =  /x5y  "  “  ^’’r'Tn  m 

direction,  i  is  the  number  of  step,  and  5y  is  the  step  size.  y„  (i  y,)  oZ).  qua 

shows  that  the  correlation  function,  p(Ay).  depends  on  the  parameter  y,  and  echo  signals  will 
be  completely  de-correlated  when  y„  £2.  As  shown  in  Figure  2,  the  function  given  by 
Equation  (7)  could  be  well  approximately  insteaded  by  a  Gaussian  function 


where  ~  2J2h^  / 

II.  experimental  procedures  and  results 

Experiments  were  performed  to  test  the  proposed  method  to  determine  ^eMive 

sgsgsislls 

Page  5 


accurate  correlation  function,  the  step  size  is  set  as  =  0.1  mm  in  our  experiment. 

At  first,  for  each  point  within  the  region  of  interest,  the  pixel  value  is  converted  to  a 

linear  scale  by  I{r)  ~  ,  where  P(F)  is  raw  B-mode  imaging  dB  data,  and  the  constant 

is  used  to  convert  from  a  0-255  log  scale  for  intensity  to  a  linear  scale.  Then,  the 
experimental  result  of  the  normalized  correlation  function  of  the  echo  signal  intensities  among 
different  image  frames  is  determined  by 


(ra  X  +  A^e,) 

= - f - 

j=i  *=i 


(8) 


here  the  average  is  over  the  region  of  interest,  e.g.  mxn  pixels.  Then,  the  correlation  function 
of  echo  signal  intensities  can  be  determined  by  using  Equation  (8).  Finally,  a  Gaussian  curve  is 
used  to  fit  the  experimental  results,  and  the  step  size,  5y,  can  be  estimated  based  on  the 
coefficient  in  that  fitting  curve. 

In  order  to  examine  the  effect  of  this  slice  positioning  technique  on  a  3-D  image 
reconstruction  and  to  further  indicate  its  accuracy,  images  of  a  tissue  immicking  phantom  from 
the  3.5  MHz  array  transducer  were  made  using  a  mechanical  positioning  system.  The  system 
was  specifically  designed  for  breast  scans,  but  has  also  been  successfully  used  to  scan  a 
transplanted  kidney  and  some  tissue  mimicking  phantoms  [9].  The  images  were  then  processed 
to  determine  the  image-based  spacing  that  would  have  been  calculated.  Post-processing  of 
images  was  implemented  in  a  modular  code  under  a  AVS  (Advanced  Visual  Systems,  Inc. 
Waltham,  MA)  and  run  on  a  Sun  Sparcstation.  Images  were  collected  from  an  ultrasound 
scanner  and  read  into  the  workstation  memory.  The  RGB  images  are  then  processed  to  select 
one  channel  for  B-mode  and  then  the  3-D  data  set  is  sliced  Md  a  single  2-D  plane  displayed  to 
select  regions  of  interest  (ROI)  to  be  processed  for  determining  the  slice  separation.  The  pixels 
contained  in  the  ROI  in  each  image  are  converted  to  a  1-D  vector  and  the  process  repeated  for 
each  image.  Each  of  these  vectors  are  then  combined  in  a  2-D  vector.  The  group  of  image 
ROIs,  which  are  at  the  same  depth  on  B-mode  images,  are  then  analyzed  to  compute  the 
correlation  between  successive  ROIs,  i.e.  ROI#I  correlation  to  ROI#2,  ROI#2  to  ROI#3,  etc. 
for  a  one  step  correlation  values.  Then  the  process  is  repeated  for  two,  three,  etc  step 
correlation.  For  this  test,  ROIs  from  groups  of  ten  images  were  used  to  determine  the  slice 
separation  for  the  center  two  slices,  e.g.,  to  position  image  6  with  respect  to  5  use  information 
from  images  0-9,  for  image  7  after  6  use  image  1-10,  etc.  This  does  assume  a  piece- wide 
smooth  motion  for  the  scan-head  but  this  is  considered  a  re^onable  first  approach  to  the 
position  estimate.  In  each  case,  the  correlation  curve  was  fit  using  a  least  squares  approach  to 
the  Gaussian  function  as  shown  in  Equation  (7').  This  information  is  then  used  to  position  the 
images  on  an  arbitrary  grid  space  and  correctly  positioned  images  can  then  be  displayed. 

A  typical  measurement  result  of  the  correlation  function  of  the  echo  intensities  as  well  as 
its  Gaussian  fitting  curve  are  shown  in  Figure  3.  Based  on  this  Gaussian  fitting  curve,  the 
optimal  value  of  the  step  size  can  be  determined  at  that  specific  depth.  Similarly,  a  group  of  the 
correlation  functions  are  given  in  Figure  4  when  the  scan-head  mechanic^ly  moved  through  the 
phantom  in  elevational  direction.  Then,  the  space  between  image  planes  is  determined  based  on 
the  coefficients  of  these  correlation  functions.  The  average  calculated  image  spacing  was 
(0. 17  ±0.01)  mm  where  the  predict  value  was  0.17  mm.  The  mean  and  stand^d  deviation  of 
the  step  size  are  calculated  based  on  several  regions  of  interest  at  the  depth.  This  indicates  that 
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the  estimating  technique  is  working  quite  well. 

III.  DISCUSSIONS 

As  shown  in  the  current  paper,  the  results  of  correlation  between  successive  images  can 
indicate  the  distance  traveled  accurately.  There  are  a  number  of  ways  to  use  this  information  to 
monitor  scanner  motion.  Many  of  these  can  be  thought  of  simply  by  reversing  the  concepts 
currently  used  to  measure  tissue  motion.  The  following  are  a  few  example  impleinenta  ions  for 
monitonng  scanner  motion  using  this  correlation  techniques  alone  or  in  the  combination  with 
Other  methods  of  motion  detection. 

1  Monitoring  tilt  and  rotation  of  transducer.  ....  -u  •  tUr. 

Measuring  different  motions  at  different  locations  within  mages  will  indicate  the 
transducer  motion  which  is  not  uniform  in  the  elevational  direction. 
transducer  slightly  can  be  detected  by  exanunmg  the  de-correlation  at 

imaee  as  long  as  tL  depth  dependent  p5/ (point  spread  function)  is  well  defined.  In  addition  to 
Sngle  traiducerSon^might  be /racked.  For  example,  several  different  ROIs  corrected 
for  the  psfin  these  regions,  could  be  monitored  such  that  differences  among  these 
nrovide^relative  motion  information.  The  difference  in  rate  among  these  images  would  indicate 
tilting  of  the  scan-head.  This  measurement  could  be  performed  across  the  entire  outlined  region 
comparing  the  value  in  the  region  proximal  to  the  transducer  to  that  of  the  corresponding  di 
value  of  the  same  lateral  location. 

tor  .motion,  i.  may  he  poyible  to 
minimize^the  computational  time  to  the  point  where  such  monitoring  can  be  performed  m  real 
time  The  advantage  would  be  that  the  position  information  could  be  computed  at  a  rate  faster 
hm'the  framrrSf  the  imaging  system.  Also  it  may  be  possible  to  set  motion  criteria  which 
will  ndicatrwSn  Icimagi  sLuld  be  taken  to  provide  an  automated  way  of  reducing  the 
*lbe?offames  acquM  the  3-D  data  set  and  fix  their  spacing  to  make  reconstruct, on 

simpler.  ,  .  ^ 

3 .  Correction  for  tissue  motion  such  as  respiratory  motion  correction  of 

These  sneckle  de-correlation  techniques  could  also  allow  for  the  correction  oi 

respiratory  motfon  when  3-D  data  set  are  reconstructed.  Elevational  as  well  as  m  plane 

correction  to  keep  a  region  of  interest  stationary  in  the  image. 

The  rate^ha/thrspeckle  pattern  is  changing  can  dso  indicate  the 
aberrators  in  the  field  of  view  for  the  scan-head.  By  looking  at  v^ous  regions  of  the  iinag  , 
sudden  or  inconsistent  regional  changes  in  correlation  rates  will  indicate  an  aberration  of  phase. 
Such  mapping  could  predict  which  data  to  throw  out  of  some  images  when  compounding  or 
reconsTructing  in  3-D  or  perhaps  when  and  where  phase  aberration  correction  would  need  to  be 

performed  in  the  images. 

‘S  Use  of  technique  with  current  coarse  positioning.  .  • 

The  currently  described  technique  could  also  be  used  for  fine  positioning  in  conjunction 
with  other  positioning  systems.  The  latter  could  monitor  the  long  range  motion  of  the  scan- 
head  but  lack  the  fine  spatial  resolution  required  for  the  best  3-D  reconstructions.  Such 
encodincT  devices  include  optical,  mechanical,  magnetic,  electric  spark  systems,  etc.  All  of 
which  have  been  used  in  the  past  for  3-D  ultrasound  imaging. 

grp“mp2r"^^^ 
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required  for  placing  the  adjacent  slice  with  respect  to  the  first.  The  de-correlation  which  then 
remained  between  the  slices  would  be  the  result  of  the  translation  in  the  elevational  direction  and 
would  be  an  improved  estimate  over  assuming  no  vertical  or  horizontal  motion. 

IV.  CONCLUSIONS 

In  this  paper,  a  method  using  the  correlation  of  echo  signal  intensities  in  B-mode 
images,  which  is  directly  related  to  the  changes  of  speckle  pattern,  in  the  elevation  direction  for 
determining  the  scan-plane  relative  positions  and  their  motions  is  described.  The  method  was 
tested  using  a  tissue  mimicking  phantom  containing  randomly  positioned  fine  solid  particles. 
The  experimental  result  was  in  reasonable  good  agreement  with  the  prediction.  However,  the 
current  method  for  the  de-correlation  of  the  speckle  can  only  be  applied  to  a  clinical  set  when 
tissue  can  be  modeled  as  a  gelatin  matrix  containing  many  small  randomly  distributed  point-like 
scatterers  with  uniform  acoustic  properties.  Nevertheless,  the  current  method  should  serve  as  a 
baseline  for  the  motion  detection  based  on  the  correlation  of  speckles. 
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APPENDIX 


We  consider  a  complex  random  process,  e.g.  the  complex  echo  signal,  including  both 
an  amplitude  and  a  phase  information.  We  are  interested  in  the  correlation  function  of  the 
complS  echo  signals  between  two  subsequent  scan  plane  positions.  The  complex  echo  signals 

are  given  by  [8] 


and 


Uir)  =  U,ir)  +  iUiir) 


U(r  +  Ar)  =  Ur{r  +  Ar)  +  it/,-  (r  +  Ar) 


where  r  and  r  +  Ar  aie  the  positions  within  a  region  of  interest  related  to  the  transducer  scan 
hpnd  The  comDletelv  random  signals  are  taken  to  have  zero  mean  and  to  have  the  same 

wtncis  of  Then  the  auto-correlation  function  of  the  resulttng 

signal  is  given  by 


p(4F)=<  U(7)U\r  +  A7)> 


(A-1) 


where  *  refers  to  the  complex  conjugation,  and  displacement  between 
planes  is  given  by  Ar .  If  this  field  is  simply  read  out  or  scanned  by  an 
point-spread  function  (psf)  Hr),  the  resulting  complex  signal  U{7)  is 

field  a{r)  by 


points  in  the  two  scan 
imaging  system  with  a 
related  to  the  original 


U{r)  =  h{r)®a{r) 


(A-2) 


where  0  is  a  convolution  operation  and  a(r)  depends  on  the  properties  of  the  medium.  The 
correlation  function  of  the  resulting  processing,  p{Ar)  is  then  directly  given  by  [7] 


p{Ar)  =  h{-Ar)  0  pa(Ar)  0  h*  (Ar) 


(A-3) 


where  p  (Ar)  =<  a{r)a  {r  +  Ar)  > .  Equation  (A-3)  shows  that  the  correlation  function 

denends  on  the  displacement  of  the  beam  related  to  scatterers  in  the  region  of  interest. 

^  A  typical  eLmple  is  the  acoustic  scattering  from  a  random  collection  ofN  scattejers  in 
and  extending  well  beyond  a  resolution  cell  or  volume.  These  N  scatterers  can  be  divided  into 

L  sets  each  set  containing  Ni  scatterers  ( I  A;  =  N),  and  all  having  the  same  displacement, 

/=! 

APi  from  the  previous  set.  In  that  case,  we  have 


N 


10,- 


(A-4) 


i=i 


and 
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(A-5) 


N 

a(r +  /lF)  =  ^|aj|e‘  ^ 

;=i 

where  a,'  is  the  magnitude  of  the  scattering  amplitude  of  the  i  th  scatterer  and  there  is  correlation 
neither  among  the  scatter  phase  in  the  first  scattering  volume  nor  among  those  in  the  second 
scattering  volume  nor  between  the  scattering  volumes.  Then  the  correlation  function  of  cjc(r)  is 

<  a(r)a*(F  +  Zir)>=5^5^<|a,'l|fly|e'^^' 

1=1  ;=1 

/=1  ^ 

since  all  phases  average  out  unless  they  match  up  identically.  Here  d{Ari)  is  the  Dirac  delta 

function.  .  .  ,-r-  j 

Finally,  the  correlation  function  of  the  readout  or  scanned  process  is  then  simplified  to 

L 

p(Ar)  =  ^pl{Ari)  (A-7) 

/=i 


where  p/  (ZlF/ )  is  given  by 


Pi  ( AF/ )  =  Ni  ||fl/  f  [h{-Ari )  ®  /i*  ( AF,  )|  .  ( A-8) 

In  the  far  field  or  the  focal  zone  of  the  transducer,  the  psf  can  be  approximately 
separated  into  the  transverse  component  due  to  the  diffraction  limited  beam,  and  the  axial 
component  due  to  the  ultrasonic  pulse 

hiAr) « (A-9) 

where  h^irj)  and  r,  are  the  transverse  components  of  the  p5/and  the  displacement  of  the 

scatterers,  and  and  are  the  axial  components  of  the  psf  and  the  displacement  of 

scatterers.  As  a  simple  example,  when  all  scatterers  are  in  motion  (relative  to  the  scan  head)  in 
the  lateral  direction,  the  correlation  function  of  the  complex  echo  signal  will  just  depend  on  the 
transducer  beam  properties  when  the  scatterers  are  in  motion  in  the  axial  direction,  that  function 
will  depend  on  the  ultrasonic  pulses. 
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abstract 


T  o  nrpvious  Daoer  fH,  a  method  for  determining  the  scan-plane 
an/o?ien.a«„^ 

speckle  on  a  series  of  B-ode  .mages  iite^ireramong  the  Images  slices 

second-order  statistics  of  echo  g  in  forming  these 

within  regions  of  i„  the  estimate  this  auto-correlation 

echo  signals  give  rise  to  an  one  J  evalnated  here  using  the  error 

Xli^a-t-ion.^X  statistical  --^aln^depends^^^^^^^^^^^^  ^^rresSmte ''th^ 

Sco"rr'ttto“fnn^  Further  elperiments  are  reported  in  the  current 


paper. 
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INTRODUCTION 

In  a  previous  paper  [1],  we  studied  a  method  to  estimate  scan  plane  motions  and  their 
relative  positions  by  using  the  de-correlation  of  the  speckles  in  ultr^onic  B-mode  images. 
Initial  results  show  that  the  method  could  be  used  to  determine  the  relahye  positions  effectively. 
The  purpose  of  this  paper  is  more  detail  to  study  the  expenmental  condiUon  which  is  required  to 
eet  more  accurate  results,  and  to  show  more  experimental  results.  ...  -  • 

^  As  given  in  the  previous  paper,  in  this  paper  it  is  assumed  that  the  statistical  properties 
of  echo  signals  follow  those  of  a  complex,  circular  Gaussian  and  the  case  is  considered  of 
diffusely  scattering  tissue  with  many  fine  particles  per  resolution  cell  [2]  and  with  no  phase 
distortion  [31.  The  method  involves  evaluating  the  second-order  statistics  of  echo  signals,  ihe 
random  processes  involved  in  forming  these  echo  signals  give  rise  an  uncertainty  in  the  estimate 
the  second-order  statistics.  This  statistical  uncertainty  is  evaluated  here  first  using  error 
propagation. 

I.  THEORY 

a.  Statistical  uncertainty  of  the  correlation  function  of  speckles 

in  B-mode  images  .  ^  u  ^  i  .  f 

The  correlation  of  the  echo  signal  intensities,  which  is  related  to  the  de-correlation  ot 

speckles  between  two  specific  positions,  e.g.  r  and  r  +  Ar  is  given  by 

R{Ar)=<I{r)Iir  +  Ar)>  (1) 

where  /(F)  and  /(F  -i-  Ar)  are  the  echo  signal  intensities  at  these  two  positions.  Ar  is  the 
displacement  between  two  points  in  the  region  of  interest.  Thus,  the  relative  scan  plane 

positions  and  their  motions  can  be  determined  from  this  correlation  function  [IJ. 

To  compute  the  auto-correlation  function  of  the  imaging  intensity  effectively  ,  we  must 
determine  how  big  of  a  region  of  interest  is  needed  to  determine  this  correlation  fonction  within 
a  specific  statistical  uncertainty.  This  includes  two  things;  1)  the  relationship  between  the 
independent  number  of  samples  and  statistical  uncertainty  of  the  auto-cortelation  fonction;  2)  the 
relationship  between  the  number  of  independent  sample  and  a  volume  of  region  of  interest. 

In  our  experiment,  we  use  the  following  equation  to  determine  the  correlation  function 

of  echo  signal  intensities 


r(AF)»4T 


(2) 


where  I{r)Iir  +  Ar)»,  Iir)»^«  I{r  +  Ar)>>,  and_«...»  stands  for 

the  experimental  sample  means  of  random  variables  /(F)/(F  +  AF)  and  /(F).  The  variance  of 
the  variable  R{Ar)  is  given  by  [4] 


(3) 


where 


(t]  =— (<  l'^{r)I^ir  +  Ar)>-<I(,r)I(r  +  Ar)>‘^) 

M 


.  2- 
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M 


-n.2. 


and 


M 

=  i_  <  {[/(r)/(F  +  Ar)-  <  /(r)/(r  +  4r)  >][/(?)-  <  I{r)  >]}  > 

M 

=  ^(2k,-r}<i>’. 

M. 

Here  M  is  the  number  of  independent  of  samples  in  the  region  of  interest  used  to  determine  the 
sample  mean  of  the  random  variables[5],  and  the  constants  k,  .  k^.  and  V  could  be  detemuned 
bv  the  relative  positions  of  two  points  r  and  r  +  4? .  Using  the  complex  Gaussian  statistma 
properties  of  the  echo  signals,  the  values  of  these  constants  can  be  determined  Md  their  values 
will  be  within  It,  e  [1,61,  il,  e  (1.31,  and  V  s  [1,2]  [6],  Equation  (A-2)  is  simplified  to 

.2 


(t|  =^(4fei -8A:2i?  +  7/?^) 


(4) 


As  an  example,  we  consider  two  limitation  cases. 

(a)  When  zlF  =  0,  we  have  =  6,  k2='^,  and  F  =  2.  Therefore,  a^-^l M, 

(b)  When  ziF  =  ~,  we  have  ki=l,  ^2  =  1’  and  F  =  1.  Therefore,  g^  =  3I  M. 

In  general  case,  we  have  0<Ar<  +°o.  Therefore,  the  variance  of  the  auto-correlation 

function  is  given  by 

(5) 

m  m 

That  means  the  variance  of  the  auto-correlation  function  is  directly  related  to  the  number 
of  indeoendent  data  samples  M.  Using  the  statistics  of  the  image  intensity,  it  is  not  difficult  to 
detemune  the  number  of  independent  data  samples  M,  which  is  related  to  the  statistics  of  the 
sample  mean  of  the  image  intensity  as  well  as  the  vanance  of  the  intensities  l/J. 

b .  Error  analysis  in  the  fitting  curve 

In  our  measurement.  Equation  (7)  in  reference  1  is  used  to  fit  the  experimental  data  p,, 

at  different  displacement  Aji  =ixy,  to  determine  the  relative  positions  of  the  scan  planes 

(called  the  step  size),  where  i  the  number  of  the  step  and  is  the  step  size.  It  is  important  to 

determine  the  uncertainties  involved  in  this  technique. 

The  difference  between  a  Gaussian  fitting  curve  and  the  expenmental  data  is  noted  by 

Uyc)  =  J,(Pi-Piof 

where  N  is  the  total  number  of  steps,  pio  is  the  measured  values  of  the  auto-correlation  function 

at  a  different  displacement  Ay^  and  p,-  =p{Ayi). 

The  optimized  step  size,  could  be  determined  from  the  condition  of  the  minimum 

L(yc).  That  is 
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(7) 


=  i(p(4y,-)-p„)^^=0  . 


/=! 


As  shown  in  the  reference  1,  a  Gaussian  function  p(Ay,)  =  is  weli  approximated  as 

the  theoretical  prediction  of  the  de-correlation  function,  where  =  2. 72ii  ■  fu  that 

case,  we  have 

(8) 


N  r 


^(}'c)  =  -2^o)'cS 
1=1 


.p.^y-vVyO. 


If  the  measured  values  of  the*  auto-correlation  function  at  different  displacements  are 
independent,  the  uncertainty  of  the  step  size  will  be  given  by 


N 

;=i 


dp. 


•{^Pjo) 


(9) 


where  8pio  is  the  uncertainty  of  the  measurement  result  of  the  auto-correlation  function  at  step 
i,  and 

dy,  [dF!  dp  Jo) 


dp  jo  (dFIdyo) 


-yJpj 


+i^Pi{Pi -Piop^^oi^yc  -l)] 

i=l  -  . 


(10) 


If  the  uncertainties  of  the  measured  auto-correlation  function  for  all  steps  are  almost  the 
same,  5p  =  5pj^  =  Pi  -  Pio,  Equation  (9)  can  be  simplified  to 

i('V) 


^  j:,.  a 


^y. 


Vy 


c  J 


\fpaoi^ylp}^i\2aoi^yl-\)p,5p] 


■5^p. 


(H) 


II.  EXPERIMENTAL  PROCEDURES  AND  RESULTS 

Experiments  were  performed  to  test  the  proposed  method  to  detemine  the  relative 
positions  of  scan  planes  by  using  the  correlation  of  the  speckles  on  B-mode  images^  A  5  U 
MHz  array  transducer,  with  the  elevation  length  is  h-6  mm  and  the  radius  curva  ur 
transducer  is  r=3.5  cm,  was  directly  used  to  collect  imaging  data  within  the  region  of  the 
interest  in  a  tissue-mimicking  phantom,  which  contains  randomly  distributed  pomt-like 
particles.  When  Ar,  >  1.7mm  the  echo  signal  will  be  completely  de-correlated,  in  order  to 

determine  the  accurate  correlation  function,  the  step  size  is  0. 1  mm  in  our  experiment. 

At  first,  for  each  point  within  the  region  of  interest,  the  pixel  value  is  then  converted  to  a 

linear  scale  [7] 

/(r) » 10'’*'^"’” 
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where  P(r)  is  raw  dB  data,  and  the  factor  F,  is  used  to  convert  from  a  0-255  log  scale  for 
intensity  to  a  linear  scale.  Then,  the  correlation  function  of  the  echo  intensities  is  given  by 


r(Ay,)^ 


(m  X  +  Ay,e,) 


;  =  1  t=l 


-1 


(12) 


Ira  n 

J=1  *=1 


Here  the  average  is  over  the  region  of  the  interest,  e.g.  mxn  pixels.  Then  Equation  (7)  of 

reference  1  is  used  to  fit  the  experimental  result  and  to  find  an  optimal  step  size  • 
reference  ‘  auto-correlation  functions  based  on  Equ^ons  (12)  and  . 

the  slice  separation.  The  image  Each  of  these  vectors  are  then  combined  in  a  2- 

dSemiie  the^slice  separation  for  the  center  two  slices  ThildoeTSme 

LmbeMsTncrea^ed  because  the  uncertainty  of  the  step  ^ae  ro'’»V'Flg“how's  t 
accuracy  of  the  mechanical  scanning  system,  is  very  small  (about  10  ^).  Fig^  J  sno«  ^ 

ScS.3ri; 
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were  recorded  at  32  fps  in  a  120  frame  cine  loop.  Results  given  in  Figure  4  asow  the 
anoearance  of  the  echogenic  cone  is  not  very  different  for  this  short  scan  when  using  either  an 
eaual  slice  separation  assuming  uniform  hand  scanning  rate  or  that  calculated  using  the 
decorrelation  rate.  However,  the  lower  graphs  illustrate  the  decorrelation  does  track  smoothly 
the  non-uniform  motion.  Black  stripes  in  the  subtraction  image  correspond  to  zeros  in  the 
integral  displacement  error  of  the  uniform  stepping  as  plotted  in  the  graph  on  the  nght.  The 
results  of  these  images  indicate  that  the  estimating  technique  is  working  quite  well  at  this  early 

stage  of  development. 

IV.  CONCLUSION: 

In  the  current  paper,  we  further  present  the  method  of  de-correlation  of  the  speckle 
pattern  in  the  elevation  direction  to  determine  the  scan  plane  relative  positions  and  their  motions. 
The  experimental  results  based  on  a  tissue-mimicking  phantom  seems  good  for  the  current 
technique.  However,  the  current  method  for  the  de-correlation  of  the  speclde  can  only  be 
applied  to  a  clinical  set  when  tissue  can  be  modeled  as  a  gelatin  matnx  contmmng  many  smal 
randomly  distributed  point-like  scatterers.  Nevertheless,  the  current  method  should  serve  as  a 
baseline  for  the  motion  detection  based  on  the  de-correlation  of  speckles. 
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Fieure  V  The  normalized  auto-correlation  function  of  echo  signal  intensity  based  on  the  region 
of  interest  in  the  focal  region  or  in  the  far-field  region.  The  solid  line  is  the  result  based  on. 
Equation  (11),  while  the  dash  line  is  the  result  of  Gaussian  fitting  curve. 

Fieure  2-  When  the  scan-plane  is  automatically  scanned  through  a  tissue  mirnicking  phantom  in 
the  elevation  direction,  the  normalized  correlation  function  of  echo  intensities  as  a  function  of 
different  number  of  steps. 

Figure  3’  When  the  scan-plane  is  manually  scanned  through  a  tissue  mimicking  phantom  in  the 
elevation  direction,  the  normalized  correlation  function  of  echo  intensities  as  a  function  of 
different  number  of  steps. 

Fieure  4-  Left  to  right  (1)  The  reconstruction  of  a  longitudinal  view  of  cone  object  in  the 
contrast  detail  phantom  from  a  3.5  MHz  array  transducer  using  assuming  equal  step  size,  (2) 
The  reconstruction  of  the  same  view  using  manually  actual  stepsize  estimated  by  the  come 
of  speckle;  (3)  The  subtract  image  between  these  two  views.  The  left  ^aph  shows  the  step  size 
changed  with  the  step  number,  and  the  right  graph  shows  the  integrated  step  size  changed 

the  step  number. 
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